Universality of the Gunn effect: self-sustained oscillations mediated by solitary waves 
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The Gunn effect consists of time-periodic oscillations of the current flowing through an external 
purely resistive circuit mediated by solitary wave dynamics of the electric field on an attached 
appropriate semiconductor. By means of a new asymptotic analysis, it is argued that Gunn-like 
behavior occurs in specific classes of model equations. As an illustration, an example related to the 
constrained Cahn-Allen equation is analyzed. 
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In semiconductors where the local current density as a 
function of the local electric field i^N-shaped, the Gunn, 
effect is an ubiquitous phenomenonEJ~El. The Gunn effectcl 
consists of time-periodic oscillations of the electric cur- 
rent flowing through an external purely resistive circuit 
attached to a semiconductor sample subject to dc volt- 
age bias. The current oscillations correspond to the gen- 
eration, one-dimensional motion and annihilation of soli- 
tary waves of the electric field inside the semiconductor. 
Besides this, the onset of the Gunn effect can be quite 
interesting, as the current may display intermittency ac- 
companied by spatio-temporal structures of the electric 
field inside the semiconductor!!!. Recently the onset of the 
Gunn instability has been analyzed by singular pertur- 
bation methods which provide the governing amplitude 
equation for long semiconductors!. Gunn-like phenom- 
ena may also explain the experimentally-observed self- 
sustained oscillations, of the current in doped weakly- 
coupled superlatticesO whose dominant transport mech- 
anisiHqis resonant tunneling between adjacent quantum 
wellstl In these cases, the oscillations are due to re- 
cycling of electric field wavefronts (charge monopolcs) 
instead of solitary wavescn The difference in the type 
of the waves may be tr.acksd to the boundary condition 
at the injecting contactliiEH. Gunn-like phenomena have 
also been numerically observed in a driiten diffusive lat- 
tice gas model of hopping conductivity^. 

A natural question that comes to mind in relation with 
these phenomena concerns their universality: Given that 
the Gunn instability appears in widely different semi- 
conductor systems and models, what are the features a 
given model has to have in order to present the Gunn in- 
stability? Notice that the Gunn effect is in principle a 
nonequilibrium phenomenon which may happen far from 
any bifurcation points. Thus the question of its universal- 
ity may not be related to linearization about fixed points 
of a renormalization transformation. Nevertheless a new 



asymptotic analysis allow us to understand deeply the 
Gunn effect and to try to give a precise meaning to the 
notion of universality far from equilibrium. This paper 
tries to give an answer to the universality question and it 
also puts the Gunn instability into perspective by com- 
paring it tor-phenomena occurring in other pattern form- 
ing systemsEJ. 

From the study of the Gunn instability in semicon- 
ductor models, we can extract the following common fea- 
tures that seem to be necessary for its occurrence: 

1 . The model should be able to support solitary waves 
moving in a privileged direction on a large enough 
spatial support. 

2. It should include an integral (over space) con- 
straint. 

3. It should have appropriate boundary conditions 
(Dirichlet, Neumann, mixed, . . . ) which render un- 
stable the stationary solutions for certain values of 
the integral constraint. 

We shall illustrate these points by constructing a sim- 
ple model that displays the Gunn instability: 
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In these equations the unknowns are u(x,t) and J(i), 
with t > and < x < L; g(u) is a function having 
a local maximum gu = g{uy[) followed by a local min- 
imum g m — g{u m ) for u > (0 < % < u m ), while K 
and (j) are non-negative parameters. Equations (|l|)-(|^) 
are to be solved with an appropriate initial condition for 
u(x, 0) > and Dirichlet boundary conditions: 
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u(0, t) = u(L, t) = P J(t), Hz. < - < M . ( 3 ) 

u m P um 

In semiconductor models u, J and 4>L correspond to the 
electric field, total current density and dc voltage bias, 
respectively. The boundary conditions (^|) correspond to 
Ohm's law relating the electric field and the current at 
the injecting and receiving contacts. (We assume that 
both contacts have identical resistivity p > for simplic- 
ity). Other boundary conditions (fixed u, mixed bound- 
ary conditions) [tlojiot qualitatively change the character 
of the solutions!! 
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FIG. 1. (a) The function 

J(t) for g(u) = 10Q{u - 0.2)(u - 0.4)(u - 0.6). Parameter 
values are K = 2, L = 100, p = 1.5 and <j) = 0.32. (b) 
The corresponding profiles of u(x, t) evaluated at the times 
marked in Part (a) of this figure. The minimum value of J 
corresponds to J* = 0, whereas the plateau at intermediate 
values of J corresponds to the solution of 2c+(J) = c_(J). 

The model represented by Equations (jl|)-(j|) with 
K = and zero-flux boundary conditions instead of 
(||) is known as the constrained Cahn-Allen equation, 
and it was recently introduced by Rubinstein and Stern- 
berg as a nonlocal reaction-diffusion model of nucleation 
akin to the mass-conserving fourth-order Cahn-Hilliard 
equationll3'E3. Equation (Q) with a fixed constant J and 
K = is the well-known bistable Fisher-Kolmogorov- 



Petrovskii-Piskunov (FKPP) equation, which includes 
among its possible solutions a variety of traveling fronts 
and pulses (solitary waves) moving on an infinite 1 D spa- 
tial supporttlrEJ. The pulses of the FKPP equation are 
unstable solutions: they either shrink or expand when 
an infinitesimal disturbance is addedO. The global inte- 
gral constraint (||) and Dirichlet boundary conditions (^) 
convert the FKPP equation in a model very similar to 
the typical semiconductor ones: the constrained Cahn- 
Allen equation. This model does not present the Gunn 
instability if K = because the x <-> —x symmetry im- 
plies no preferred direction of motion for traveling waves. 
A large enough nonzero convective term K > breaks 
the x <-> — x symmetry and it privileges waves moving 
from left to right. The resulting model satisfies the con- 
ditions 1 to 3 above and it displays the Gunn effect; see 
Fig. 0, It may be observed that the present model is 
also related to Kroemer's model of the Gunn effect in 
n-GaAsa: we just change the convection coefficient to a 
constant K in Ampere's law and set the diffusivity equal 
to one in the dimensionless Kroemer's model studied in 
Ref. [l2| These changes exclude the straightforward ex- 
tension of our previous asymptotic analysis, as we cannot 
use the shock waves and particular solutions-specific of 
Kroemer's model to describe the Gunn effeciJla. 

To understand these results, we shall assume that 
e = 1/L <C 1. Then it is convenient to rewrite Equa- 
tions (|l|)-(||) in terms of the 'slow' variables s — et and 
y = ex. The result is 
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In the limit e — > the solutions of this system are 
piecewise constant: on most of the y-interval u is equal 
to one or another of the zeros of g(u) — J, separated 
by transition layers that connect them. At y = and 
y = 1 there are boundary layers (quasi-stationary most 
of the time), which we will call injecting and receiving 
layers, respectively. Let us assume that uu < 4> < u m 
and denote by U\{J) < u 2 (J) < us(J) the three zeros 
of g(u) — J. Let the initial profile u(y,0) satisfying (|^) 
be a square bump u — uz{J) for Yi(0) < y < 5^(0) 
and u = u\(J) elsewhere plus terms of order e, as in the 
time marked by (1) in Fig. 0(b). Located at y = Yi and 
y = Y2, Y\ < Y2, there are sharp wavefronts of width 
O(e) connecting u = u\{J) and u — u^{J). This initial 
profile will naturally evolve into the Gunn effect as time 
oes on (see below). The initial value of J follows from 



= Ul (J) + [u 3 {J) - ui(J)] (Y 2 - Y^ + O(e). (6) 

The boundary layers and the fronts connecting u\(J) 
and Us( J) are built from trajectories of the phase plane: 



2 



du 



dv 



+ g{u) - J, 



(7) 



where £ = e~ l [y — Yi(s)], c = dYi/ds and p, = K — c. 
The boundary layers are separatrices connecting the ver- 
tical line u — pJ in the phase plane (u, v) to the saddles 
(tti,0) or (m3,0) for c = 0: — > ttj(J) as a; — * oo 
and u(a;) — > u,-(J) as (x — L) — » — oo (i = 1,3) are the 
matching conditions. For each fixed value of J between 
g m and gu we can find a unique value c+(J) such that 
u{— oo) = u\(J) and u(oo) = i«3(J) [corresponding to a 
heteroclinic orbit connecting (ui, 0) to (113, 0) with v > 0] 
and a unique value c_(J) such that u(— 00) = u^{J) and 
u(oo) = u\(J) [a heteroclinic orbit connecting (113, 0) to 
(tiljO) with v < 0]. The functions c±(J) are depicted in 
Fig. 0. They intersect when J = J* given by 



J* = 



1 



u 3 - u x 



g{u) du, c± = K. 
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FIG. 2. The functions c±(J) for K = 2, and p(u) as in Fig. 
1. When J c > (vertical line at J = J T ) only one wave 
is shed during each period, whereas multiple shedding may 
occur for J* < J c < J . 

Starting at s — 0, the fronts li(s) move with speeds 



ds 



= c+(J), 



dY 2 
ds 



c_(J), 
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whereas their positions are related to the bias <f> through 
(^). We find an equation for J by differentiating (@) and 
then inserting (ph in the result: 



dJ 
ds 



= A(J) [ C+ (J)-c_(J)], 
A= } U3 - Ul) \ >0, 



(10) 
(11) 



where g[ = g'(ui) and we have used that g(ui(J)) = J 
implies dui/dJ = \/g'{ui). This is a simple equation 
for J demonstrating that J tends to J* exponentially 
fast. Notice that this is a very simple explanation of the 
well-known observation that a pulse detached from the 
boundaries moves at constant speed and J, given by the 
equal area rule (||).ld 

After a certain time, the wavefront Yi reaches 1 and 
we have a new stage governed by (0) with Y-x = 1 and 
Y\ given by ([)]). The equation for J becomes dJ/ds = 
Ac + > and its solution increases [compare J and it at 
time (2) in Fig. [ij until it surpasses the value J c such 
that U2(J) = pJ- (At J c , [du/dx] x =o changes sign aiid 
the quasistationary injecting layer becomes unstable )Ej. 
Let si be the earliest time at which J = J c . After 
s = s\, the profile of u changes within the boundary 
layer at y = 0: this injecting layer becomes unstable and 
it sheds a new wave during a fast stage described by the 
time scale r = (s — si)/e. To find what happens next 
we need to perform a more complicated analysis keep- 
ing 0(e) terms in the outer (bulk) expansion of u and J 
and just the leading-order term in all inner expansions 
(boundary layers and wavefronts). This calculation-has 
been performed in detail for a semiconductor modelliS It 
can be shown that the shedding of a new wave from the 
injecting layer is governed by the following semi-infinite 
problem for x > 0, —00 < r < 00: u(x, r) (far from the 
old wave dying at y = 1) solves (Q) and w(0, r) = pJ(r; e), 
with J(r;e) = J c + £JW(t), 



j( 1 )(r) = /i'(r)+a/i(T)- 7 / e~^ T -^ h(t) dt, (12) 



ft(r) = (it 3 - tii) c+ (r - r Q ) - / [u(x, r) - u x ) dx, (13) 

Jo 



(in this equation all functions of J are calculated at 
J = J c ; To pis a constant and a, (3 and 7 are positive 
parameters)lia and the following matching condition on 
an appropriate overlap domain: u(x, t)— uo(x; J(s)) <C 1, 
as t — > —00, s — > s±—. Here Mo(a;; >/(s)) is the quasista- 
tionary injecting layer solution of (j?)) with p = K such 
that ito(0; J(s)) = pJ{s) and tto(oo; J(s)) = Ux(J(s)) for 
s < si, J(si) — J c - The function h{r) is the area lost due 
to the motion of the old front during the time t minus 
the instantaneous excess area under the injecting layer. 

The solution of the previous semi-infinite problem re- 
veals the formation, growth and motion of a new pulse 
in the injecting layer, driven by h(r) through the effec- 
tive excess current (Q). This process ends when the 
new pulse is bounded by two well-formed wavefronts (de- 
tached from the injecting layer) which are located at I3 
and Yi, Y3 < Yi [see the u-profile at time (3) in Fig. 
|l|(b) in which Y3 and Y4 have already moved from their 
initial positions O(elne) at the beginning of this stage]. 
It may be seen that the injecting layer becomes unstable 
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and sheds a new wave when its width reaches a critical 
size Ay = 0(elne)H. 

If 4> is large enough, we have a stage where the old 
wavefront located at Y\ < 1 coexists with the newly 
formed pulse bounded by the two wavefronts located at 
y 3 and Y 4 : 

<t> = ui(.7) + MJ) - ui(J)] {l-Yx+Yi- Y 3 ) 

+0(e). (14) 
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FIG. 3. (a) Density plot for u(x,t) with p = 4, <j> = 0.3 and 
gr(u) as in Fig. 1 (lighter color means larger u). Here multiple 
shedding of pulses occurs at the injecting layer: two pulses 
are formed during each period. The secondly shed reaches 
and overtakes the first one. (b) The corresponding profile of 
J it). 

Differentiating this equation and using that Y\ and 
Y3 move with speed c + whereas Y4 moves with speed 
c_, we obtain dJ/ds = A (2c + — c_). Starting from 
J c , J decreases further to JT [the zero of (2c+ — C -)L if 
2c + (J c ) < c_(J c ) (the stable case with J c > jt in Fig. g). 
After the old wave reaches y = 1, we get again Equations 
@-© and recover the initial situation. Thus a full pe- 
riod of the Gunn oscillation is described; see Fig. [|. On 
the other hand, if 2c+(J c ) > c_(J c ) (J* < J c < J 1 *), 
J increases after the formation of the new pulse and 



it is possible for the injecting layer to shed more waves 
into the bulk as shown by the numerical simulations of 
Fig. 0. How many waves are shed depends both on the 
value of J c (and therefore on the injecting resistivity p) 
and on the length L. A rough estimation would give 
(n + l)c + (J c ) = rtc_(J c ) as e — > for the number n 
of shed waves. This shedding mechanism seems to have 
the effect of breaking the spatial coherence of the sample 
which may lead to complex spatio-temporal phenomena 
(intermittencies with a varying number of pulses present 
in the sample at different times). The unstable case will 
be further analyzed in the near future. 

In conclusion we have investigated what are the main 
features that a given model should have in order to 
present the Gunn effect. These features are demonstrated 
by studying a simple model by means of a general asymp- 
totic analysis corroborated by direct numerical simula- 
tions. As a result the Gunn effect is reduced to solving 
a sequence of very simple problems (one equation for J 
each time) plus a canonical problem for shedding new 
pulses. Our asymptotic analysis explains qualitative and 
quantitativally the formation, motion and annihilation 
of pulses in the Gunn effect. This work sheds light on 
several puzzling aspects of the Gunn oscillations (see the 
Chapter on open problems in Ref. |2^) : (i) why do pulses 
move with the well-known equal-area- rule velocity at con- 
stant J when they are far from the contacts? [the cor- 
responding current is a stable equilibrium of (fffl]; (ii) 
how does the wavespeed change when it arrives to the 
receiving contact?; (iii) how are new waves created at 
the injecting contact? In addition we have described a 
new instability mechanism consisting of multiple pulse 
shedding during each oscillation of J, which appears for 
appropriate values of the boundary parameters at the 
injecting contact. Similar work has been performed in 
diverse semiconductor models: Gvwwi oscillations in ul- 
trapurc closely compensated p-GetaLKroemer's model 
of Gunn oscillations in bulk n-GaAal3 and slow oscilla- 
tions in semi-insulating GaAsa. A modification of the 
asymptotic method presented here describes the charge 
monopole recycling responsible for the_self-oscillations in 
n-doped weakly coupled superlatticesEU. Irrespective of 
the physical mechanism responsible for the existence of 
wavefront and pulses, our asymptotic method describes 
the Gunn oscillations in these models. The model pre- 
sented here perhaps illustrates in the simplest way what 
the method consists of: (i) find the equations and bound- 
ary conditions which characterize the shape of the wave- 
fronts and their speed as functions of the current density 
J; (ii) derive the equations which determine J as a func- 
tion of the slow time scale depending on the number of 
wavefronts present in the sample. The field profile fol- 
lows adiabatically the evolution of J; and (iii) add the 
semiinfinite problems responsible for wave shedding at 
the contacts. Solution and matching of these problems 
yields an approximation to the Gunn effect in the given 
model. Of course solving some of these steps may be in 
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itself a rather complicated technical problem for partic- 
ular models requiring special asymptoticsEj. 
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